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ABSTRACT 

Energetic feedback from supernovae (SNe) and from active galactic nuclei (AGN) are 
both important processes that are thought to control how much gas is able to condense 
into galaxies and form stars. We show that although both AGN and SNe suppress star 
formation, they mutually weaken one another's effect by up to an order of magnitude in 
haloes in the mass range for which both feedback processes are efficient (10 1125 M Q < 
wi200 < 1O 125 M0). These results demonstrate the importance of the simultaneous, 
non-independent inclusion of these two processes in models of galaxy formation to 
estimate the total feedback strength. These results are of particular relevance to semi- 
analytic models, which implicitly assume the effects of the two feedback processes to 
be independent, and also to hydrodynamical simulations that model only one of the 
feedback processes. 
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1 INTRODUCTION 

Feedback from the formation of massive stars, probably in 
the form of large-scale winds driven by supernova (SN) ex- 
plosions, and from accreting black holes (BHs) associated 
with active galactic nuclei (AGN), are thought to regulate 
the star formation rates (SFRs) of low - and high-ma ss galax- 
ies, respectively (for a review see e.g. lBensonll2010l ). Obser- 
vationally the two feedback effects are frequently seen to 
be coeval. In a global sense, the evolution of the cosmic 
SFR and the luminosity densit y of quasars are tightly corre- 
lated (|Bovle fc Terle vich 1998) and on the scale of individual 
objects post-sta rburst galaxies are fou nd to preferentially 
host active BHs (|Kauffmann et al.ll2003l ). Both AGN and SN 
feedback are now routinely included in semi-analytic and nu- 
merical simulations. However, it is freq uently assumed, par- 



ticularly in semi-ana l ytic models (e.g. iBenson et al. 2003; 
|Pe Lucia et"afl |2004 bower et all l200d : iGuo et all I201J) . 



that the effect of one form of feedback is unperturbed by 
the inclusion of the other. 

In semi-analytic models it is generally assumed that the 
amount of gas reheated by SNe is oc rh,v a , where m* is the 
galaxy star formation rate, v is some characteristic velocity 
associated with the galaxy and a is a free parameter. BHs 
grow primarily in major mergers, but there is no feedback 
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associated with this 'quasar' mode, as its energy is just as- 
sumed to be part of the SN feedback energy that is released 
by the associated starburst. AGN feedback is included in 
these models by coupling the BH directly to the halo with- 
out regard for the galaxy, either by assuming that the BH 
growth decreases the cooling rate in the hot halo in pro- 
portion to the BH accretion rate l|Guo et al.ll201ll ). or by 
shutting off cooling if the BH is massive enough for the Ed- 
dington luminosity to exceed some mu ltiple of the cooling 
luminosity of the halo (jBower et al.ll2006l) . 

It is by no means a trivial assumption that SN and AGN 
feedback act independently of each other. Both feedback 
processes redistribute gas inside of galaxies in a complex 
and non-linear way. It is possible that the effect of one feed- 
back process, e.g. the factor by which the cooling of hot halo 
gas is reduced or the amount of cold gas that is heated or 
ejected, depends on the pr e sence of ot her feedback processes . 
Indeed, iPawlik fc Schavd (|2009l ) and iFinlator etafl (|201ll ) 
demonstrated that SN feedback and photo-heating, which 
semi-analytic models also assume to act independently, am- 
plify each other's suppression of the SFR of low-mass galax- 
ies. 

The aim of this letter is to use self-consistent hydrody- 
namical simulations to investigate the mutual amplification 
of AGN and SN feedback. We will show that in haloes in 
the mass range 10 11 ' 25 M < m 20 o < 10 125 M Q AGN and 
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SN feedback act to suppress one another's effects by up to 
an order of magnitude. 

This letter is structured as follows: In Sec.[2]we describe 
the numerical methods employed in this work. In Sec. [3] we 
quantify the effects of our models for SN and AGN feedback 
on the SFR in haloes of different masses and in Sec. [4] we 
discuss our findings and conclude. 



2 METHOD 

This study is based on smoothed particle hydrodynamics 
(SPH) simulations of representative volumes of the Uni- 
verse. Gravitational forces and the equations of hydrody- 
namics are solved using a significantly extended version of 
the parallel PM Tree-SPH code GADGET ill (last described in 
ISpringell 120051 ). a Lagrangian code used to calculate forces 
on a particle by particle basis. All of our simulations assume 
a ACDM cosmology with parameters determined from the 
3-yr Wilkinson Microwave Anisotropy Probe (WMAP) re- 
sults, fi m = 0.238, Ov = 0.762, D. h = 0.0418, h = 0.73, 
as = 0.74 and n B = 0.9 51. These values are c onsistent with 
the 7-yr WMAP data JKomatsu et alj|201lh . 

In addition to treating gravitational and hydrodynamic 
forces, the simulations need to follow the galaxy formation 
processes that happen on small scales. The simulations track 
star formation, SN feedback, BH growth and AGN feed- 
ba ck, radiative cooling a nd chcm odynamics, as described 



Schaye fc Dalla Vecchial d2008h . iDalla Vecchia fc Schavel 
ll2008h . lBooth fc Schavel lJ2009h . lWiersma et all (|2009al ) and 



respectively. This physical model 
the OWLS suite of simulations 



IWiersma et al.l ((2009b) 
is denoted 'AGN ' in 
l|Schave et al.ll2010t ). 

For the purposes of this work, our prescriptions for SN 
and AGN feedback are the most important aspects of the 
physical model so we describe these in some detail here. 
Feedback from SNe is implemented by injecting approxi- 
mately 40% of the energy released by Type II SNe locally as 
kinetic energy. The rest of the energy is assumed to be lost 
radiatively. Each newly formed star particle kicks on aver- 
age 2 of its neighbouring gas particles into the wind. The 
initial wind velocity is 600 km/s, w hich is consistent w ith 
observations of galaxy outflows (e.g. IVeilleux et al.ll2005l ). 

BH growt h and AGN feedback i s implemented using 
the method of iBooth fc Schavel (120091) w h ich is , in turn, a 
modification of that from ISpringel et al.l (|2005r ) . We regu- 
larly run a friends-of-friend halo finder and insert seed mass 
BHs (m sccd = 10 5 M Q ) into every halo of mass > 10 10 M 
that does not yet contain a BH. These seed BHs then grow 
both through merging and gas accretion (which is limited 
to the Eddington rate). Accretion rates in low-density gas 
(rtH < 10~ x cm" 3 ) are assumed to be equal to the Bondi- 
Hoyle rate. For higher- density, star-forming gas the Bondi- 
Hoyle rate is boosted by a factor (tih/10 -1 cm 3 ) 2 to com- 
pensate for the lack of a cold, interstellar gas phase and the 
finite resolution (see Booth & Schaye 2009 for a full discus- 
sion). The BH growth rate is related to its accretion rate 
by jtibh = (1 — e r )m accl ., where e r = 0.1 is the radiative 
efficiency of the BH. 

1 The parameter ag is 2<r lower in the WMAP 3-year data than 
allowed by the WMAP 7-year data. 



The amount of energy available for AGN feedback is 
then given by E = efe r m accr c 2 , where c is the speed of light 
and £f is a free parameter, the 'feedback efficiency'. The 
BH builds up a reservoir of energy until it is capable of 
heating one gas particle by a temperature AT = 10 8 K. 
This energy is then injected thermally into one of the BH's 
neighbours, chosen at random. The feedback efficiency is 
tuned to reproduce the global BH density at z = 0. In the 
fiducial runs et = 0.15. 

The fiducial simulations analysed in this work are per- 
formed in cubic volumes of 50 co-moving Mpc/h and con- 
tain 256 3 particles of both gas and dark matter (DM). For 
the assumed cosmological parameters this corresponds to 
DM and (initial) gas particle masses of 4.1 x 10 8 M.®/h and 
8.7 x 10 7 Mq//i, res pectively. Initial conditions are generated 
with CMBFAST (jSeliak fc Zaldarriaga|[l996T ). and evolved 
linearly to the simulation starting redshift of z — 127. The 
comoving gravitational force softening is set to 1/25 of the 
initial mean interparticle spacing and is limited to a maxi- 
mum physical scale of 2 kpc/h, which is reached at z ~ 3. 

Each simulation is run four times, starting from the 
same initial conditions. One realisation includes no feedback 
processes, the next two include either AGN or SN feedback. 
The final run includes both AGN and SN feedback. Simu- 
lations with this last physical model and the same resolu- 
tion as used here reproduce the observed z = relations 
between BHs and the mass and velo city dispersion of their 
host galaxies l|Booth fc Schavell2009l ) as well as the observed 
optic al and X-ray propertie s of the groups in which they re- 
side (JMcCarthv et alj|2010h . 

In the simulation where we treat AGN feedback but 
not SN feedback, the AGN grow over-massi ve when com- 
pared to the observed global density of BHs (jShankar et al.l 
1 20041) and the z = relations betwe en BH mass and galaxy 
bulge mass (e.g. Haring fc Rixll2004l) and velocity dispersion 



(Tremainc et al 



2002! ) . We therefore run one additional sim- 



ulation in which the AGN efficiency is increased by a factor 
6.87 (the ratio of the global BH densities in the simulations 
with and without SN feedback) to reproduce the scaling re- 
lations even in the absence of SN feedback. 

For our analysis we use snapshots of the simulations, 
which are saved at discrete output redshifts with interval 
Az = 0.125 at < z < 0.5, Az = 0.25 at 0.5 < z < 4 and 
Az = 0.5 at 4 < z ^ 6. At each of these outputs we use a 
Friends-of-Friends halo finder with linking length b = 0.2 to 
obtain a list of DM haloes. Their masses are then determined 
using the SubFind code (|Dolag et al. 2009), which employs 
a spherical-overdensity criterion centered on the most bound 
particle in each halcQ We limit our analysis to haloes of mass 
m.200 > 10 10 - 75 M , corresponding to > 100 DM particles. 



3 RESULTS 

We begin by considering the effect of each feedback process 
on the cosmic SFR. Each panel in the top row of Fig. [T] 
shows the contribution of haloes in a given mass range to 



2 All halo masses quoted in this work are spherical overdensity 
masses defined as the mass within a sphere that encloses a mean 
density 200 times the critical density of the Universe. 
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Figure 1. Effect of SN and AGN feedback on the evolution of the cosmic SFR density. Different colours represent different feedback 
models. Purple curves show simulations that include no feedback; red curves include both AGN and SN feedback and orange (blue) 
curves show simulations that include only SN (AGN) feedback. Top panels: the contribution to the cosmic SFR density from haloes 
with different masses as a function of redshift. Bottom panels: the factor by which the SFR is decreased relative to the simulation that 
includes no feedback. The horizontal, grey line in each of these panels shows no-suppression. The blue, dotted lines demonstrate that 
changing the AGN feedback efficiency has almost no effect on the overall SFRs. At the low-mass end (m2oo < 10 1125 Mq) SN feedback 
dominates and AGN feedback dominates at the high-mass end (rri200 > 1O 12 M0), but at intermediate masses both feedback processes 
contribute to the decrease in the SFR. 



the cosmic SFR as a function of redshift in each of the simu- 
lations. The colour denotes the physical model: no feedback 
(purple), SN feedback only (orange), AGN feedback only 
(blue) and both AGN and SN feedback (red). We denote the 
global SFR in the simulation that includes no feedback with 
p, , and the simulations that include only SN feedback, only 
AGN feedback and both forms of feedback as p*,sn, P*,agn 
and p*,sn+agn, respectively. We can then define suppression 
factors, Sx = -^— , where X represents one of the subscripts 
introduced above. The bottom row of Fig. [TJ shows the sup- 
pression factors corresponding to each of the upper panels. 
In each of the lower panels, the horizontal, grey line shows 
no-suppression (i.e. Sx — 1). 

In the least massive haloes (m2oo < 10 11,25 Mq; left 
panels) AGN feedback has not yet become efficient so the 
simulations with no feedback (purple) and AGN feedback 
(blue) lie very close to one another. The simulations with 
SN feedback (orange) and both AGN and SN feedback (red) 
also nearly overlay one another. This indicates that in this 
mass range, SN feedback dominates the suppression of the 



SFR. Conversely, at the highest masses (171200 > 10 12 Mq; 
right panels) SN feedback becomes increasingly inefficient 
and now the simulations that include no feedback (purple, 
solid curve) and only SN feedback (orange) lie close to one 
another, while the AGN only (blue) and AGN and SN feed- 
back (red) lines also lie very close to one another, indicating 
that in this mass range SN feedback does little to suppress 
the SFR. At intermediate masses (10 1125 M Q < m 2 oo < 
10 12 M© ; middle panels) the two feedback processes suppress 
the SFR by comparable amounts, although at low redshift 
AGN feedback becomes less efficient. 

In the absence of SN feedback, the BH population grows 
to be much more massive than observed. We therefore con- 
sider one extra simulation in which et is increased by a factor 
of 6.87 (the ratio of the global BH densities in the simula- 
tions with and without SN feedback) in order to bring the 
total mass in BHs back in line with observations. The dotted, 
blue curves in the bottom panels of Fig. Q] show suppression 
ratios in this simulation. It is notable that this line is almost 
identical to the AGN only simulation with the standard ef- 
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ficiency, indicating that in both cases the AGN injects the 
same amount of energy into its surroundings and suppresses 
the SFR by the same factor. This occurs because if the BH 
feedback efficiency is increased by some factor, the total 
mass in BHs decreases by nearly the same factor, so that 
the total amount of energy output is almost independent of 
ef. This suggests that BHs are growing until they have out- 
put some critical amount of energy, at whi ch point they are 
capable of r egulating their own growth fsee lBooth fc Schavd 
120091 . I2OIC1 for a full discussion of this point). We therefore 
do not show the high-efficiency simulation in the rest of this 
letter but note that our results are insensitive to this choice 
and that the only effect of the feedback efficiency is to scale 
the BH masses. 

In the top panel of Fig. [2] we show the 2 = suppression 
factors as a function of halo mass for SN feedback (orange) , 
AGN feedback (blue) and both feedback mechanisms (red). 
SN feedback accounts for most of the suppression of the SFR 
in haloes with masses ^ 10 12 Mg, but its efficiency falls to 
almost zero for larger halo masses. Above 111200 — 10 12 Mg, 
AGN feedback accounts for the majority of the suppression 
of the SFR. We have verified that these same trends hold at 
higher redshift. 

In order to quantify the mutual amplification of SN 
and AGN feedback, w e define an 'amplification factor' 
IjPawlik fc Schavel2009r ) for the two feedback processes using 
their suppression ratios 



100 



SsN+AGN 

Ssn x Sagn 



(1) 



A value x — 1 indicates that AGN and SN feedback suppress 
the SFR independently of one other. A value x > 1 (X < 1) 
indicates that AGN and SN feedback amplify (weaken) each 
other's ability to suppress the SFR relative to the case where 
they act independently. 

The bottom panel of Fig.[2]shows the amplification fac- 
tor, \, as a function of halo mass. At both low (771200 < 
10 11 M Q ) and high (m 2 oo > 10 12 ' 5 M Q ) halo masses, the 
mass ranges where only one of the two feedback processes 
are effective, the amplification factor x ~ 1; indicating that 
each process operates independently of the other. However, 
in the intermediate mass range j(<l, indicating that the 
effect of including both AGN and SN feedback is almost an 
order of magnitude weaker than it would be if they each re- 
duced the SFR by the same factor as when they act in iso- 
lation. 

Eq. [T] holds true if the factor by which each feedback 
process suppresses the SFR is independent. If, however, the 
feedback effects are additive and the amounts by which they 
decrease the SFR are independent then the equation takes 
on a slightly different form. The conclusion that the two 
feedback processes weaken each other also holds if they com- 
bine additively instead of mutliplicatively. This can easily be 
seen by considering that if SNe and AGN each reduce the 
SFR by more than a factor of 2, which is the case for much 
of the redshift and mass range, then their combined effect 
would yield a negative SFR if the two feedback processes 
were independent and combined additively. 

We have verified that decreasing the simulation volume 
by a factor of eight, while keeping the resolution unchanged, 
has a negligible effect on the results. If we decrease the 
numerical resolution then the magnitude of the feedback 
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Figure 2. Top panel: Present-day star formation suppression fac- 
tors as a function of halo mass for the simulations that include 
SN feedback (orange), AGN feedback (blue) and both types of 
feedback (red). Bottom panel: Amplification factor, \ (Eq. [TJ, 
as a function of halo mass. Values of x = 1 indicate that SN 
and AGN feedback operate independently of one another. Val- 
ues x < 1 indicate that SN and AGN feedback mutually suppress 
one another's effects. In haloes with masses I0 11 — I0 12 ' 5 Mg, the 
two feedback effects weaken one another's impact on the SFR by 
almost an order of magnitude relative to the case in which they 
operate independently. 



suppression becomes smaller, although the same qualita- 
tive trends hold as for the fiducial case, with SNe (AGN) 
dominating the suppression in low- (high-) mass objects, 
and an intermediate-mass regime where the two processes 
weaken one another's effects. The halo mass at which AGN 
feedback begins to effectively suppress the SFR is higher 
for the low-resolution simulation because in this simula- 
tion gas densities at the centres of haloes are significantly 
lower than in the high-resolution case, so BHs grow more 
slowly and begin to a ffect galaxies at a later time (see also 
iBooth fc Schavd 120091 ). Thus, increasing the numerical res- 
olution only strengthens our main conclusions. 



4 DISCUSSION 

Energetic feedback from SNe and from AGN are both im- 
portant processes that are thought to control how much gas 
is able to condense into galaxies and form stars. Using a 
set of cosmological SPH simulations that include, amongst 
other ingredients, both of these feedback processes, we have 
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investigated how each (and both) of these processes affect 
the star formation history of the Universe. We demonstrated 
that, in the models, both AGN and SNe suppress the cosmic 
SFR when modelled in isolation, with SNe suppressing SF 
primarily in low- mass objects (771200 < 10 12 Mq) and AGN 
feedback operating primarily in massive (m.200 > 10 12 Mq) 
objects. However, if both AGN and SN feedback are included 
in a simulation, the factor by which the SFR is suppressed 
is significantly smaller than would be expected if the two 
processes had an independent effect on the SFR. This oc- 
curs in the regime where both AGN and SN feedback are 
able to suppress the galaxy SFR coevally, i.e. in haloes of 
intermediate mass (10 1125 M Q < m 2 oo < 10 12 M ). The net 
effect is that they weaken one another's ability to suppress 
the galaxy SFR. 

We caution the reader that our simulations have not 
fully converged with respect to resolution and that the halo 
mass at which BHs begin to affect the galaxy SFR changes 
with resolution. Our simulations contain neither the numer- 
ical resolution or the physics to model the multi-phase in- 
terstellar medium. The factor by which AGN and SN feed- 
back suppress the SFR is therefore not well converged. How- 
ever, we verified that increasing the numerical resolution 
only strengthens the effects we discuss here. We note also 
that there is significant freedom in how we choose to model 
feedback processes in our simulations, and that taking sig- 
nificantly different choices about how and when energy is 
injected into galaxies may affect the results obtained. 

However, our finding that AGN and SN feedback 
weaken each other's effect on the SFR should not come as 
a surprise. Both SF and BH growth are thought to be self- 
regulated: feedback injects sufficient energy for the outflow 
to balance the accretion rate when averaged over sufficientl y 
large length and time scales (e.g. iBooth fe Schaveil2010f ). 
Since BH accretion episodes are generally accompanied by 
nuclear SF, both feedback processes will contribute to the 
nuclear outflow. If SN feedback is turned off, then the BH 
will compensate by injecting more energy. This extra AGN 
feedback will not only limit the BH's own growth, it will 
also reduce the SFR. Hence, we expect the reduction of the 
SFR due to AGN feedback to be greater in the absence of 
SN feedback. 

Our work complements that of lPawlik fc Schavel f2009), 
who found using similar methods that SN feedback and 
photo-heating by the UV background amplify one another's 
effects on the SFR in low-mass galaxies. SN feedback 'puffs 
up' galaxies, making it easier for them to be evaporated 
by the UV background, and in turn the UV background re- 
moves gas from the outskirts of the galaxy, making it easier 
for SN feedback to drive out more gas. The general conclu- 
sion that we can draw from these results, is that the way 
in which any one feedback process redistributes gas around 
the galaxy has the potential to either amplify or suppress 
the ability of other feedback processes to suppress the SFR. 
It is therefore vital that studies treat all feedback processes 
simultaneously and in a non-independent manner. 

These results have important implications for semi- 
analytic models, which generally make the implicit assump- 
tion that all types of feedback act independently of one an- 
other. Additionally, hydrodynamic simulation s that model 
AGN in the absence of SN-driven winds (e.g. iFabian et al.l 



[2010J) may draw incorrect conclusions regarding the effect of 
AGN feedback 
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